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We study the spin- 1/2 Heisenberg model in a lattice that interpolates between the triangular and 
the kagome lattices. The exchange interaction along the bonds of the kagome lattice is J, and the 
one along the bonds connecting kagome and non- kagome sites is J', so that J' = J corresponds to 
the triangular limit and J' = to the kagome one. We use variational and exact diagonalization 
techniques. We analyze the behavior of the order parameter for the antiferromagnetic phase of 
the triangular lattice, the spin gap, and the structure of the spin excitations as functions of J'/ J. 
Our results indicate that the antiferromagnetic order is not affected by the reduction of J' down 
to J'/ J ~ 0.2. Below this value, antiferromagnetic correlations grow weaker, a description of the 
ground state in terms of a Neel phase renormalized by quantum fluctuations becomes inadequate, 
and the finite-size spectra develop features that are not compatible with antiferromagnetic ordering. 
However, this phase does not appear to be connected to the kagome phase as well, as the low-energy 
spectra do not evolve with continuity for J' — > to the kagome limit. In particular, for any non- 
zero value of J', the latter interaction sets the energy scale for the low- lying spin excitations, and a 
gapless triplet spectrum, destabilizing the kagome phase, is expected. 

PACS numbers: 75.10.Jm, 75.40.Mg, 75.50.Ee 



INTRODUCTION 

Geometrically frustrated antiferromagnets are the his- 
torical candidates for the realization of a spin liquid 
ground state. Indeed, the spin-half Heisenberg antifer- 
romagnet on the triangular lattice was the first model 
to be proposed by Anderson and Fazekas 0,0 m 1973 
as a system where geometric frustration and quantum 
fluctuations could prevent zero-temperature magnetic or- 
dering in two dimensions, stabilizing instead a ground 
state with gapped spin excitations and exponentially 
decaying correlations. Since then, a good amount of 
work has been devoted to investigate the nature of 
the ground state of the triangular Heisenberg model 
H MM HHSEillllHHmH which remained an 
open question until quite recently. At present, however, 
there is a general consensus on the existence of long-range 
antiferromagnetic order following a 120° Neel pattern in 
the ground state of this model ja, 0, El 0] : frustration 
and quantum fluctuations on the two-dimensional trian- 
gular lattice are not strong enough to stabilize a non- 
magnetic ground state. 

A more promising candidate for a disordered ground 
state can be obtained through a 'dilution' of the triangu- 
lar lattice, leading to the so-called kagome net (Fig. QJ. 
In fact, on this geometry, due to the lower coordination 
(z = 4 compared with z — 6 in the triangular case), frus- 
tration is much stronger and even in the classical limit it 
gives rise to an infinite number of classical ground states, 
with ordered and disordered configurations degenerate 
in energy. [lH UtI Hs | Due to the extensive entropy of the 



classical ground state, the so-called order from disorder 
mechanism - usually stabilizing, among degenerate man- 
ifolds, long-range ordered configurations - is much less 
effective than in other frustrated models. In particular, 
while harmonic fluctuations select planar configurations, 
they turn out to be completely insensitive to their de- 
gree of order, and only non-linear effects eventually sta- 
bilize a classical ground state with Neel correlations with 
a -\/3 x \/3 pattern (see Fig. 1). Whether such classical 
minimum energy configuration possesses true long-range 
order h& or it is a critical point with power-law correla- 
tions 0, is still an open issue. In any case, due to 
the particularly delicate mechanism leading to the classi- 
cal antiferromagnetic order, the latter is expected to be 
easily destabilized by quantum fluctuations. 

The investigation of the quantum Heisenberg antiferro- 

! has been capturing increas- 

12! El El EH E3, El El 

for some time now. Mainly on the basis of numerical 
work, evidence has accumulated supporting a spin-liquid 
ground state, even though a very peculiar one: this, in 
fact, would be characterized by a small gap (~ J/20) to 
spin excitations, and by an exponential ly large number of 
singlets contiguous to the ground state. |M I27I El lil l3f| 
The classification of spin liquid of Type II has been re- 
cently proposed [2!j to classify this particular behavior. 

On the experimental side, the triangular geometry 
provides the scenario for interesting physical phenom- 
ena taking place in the superconducting compounds 



Na x Co02- v H 2 



and the organic materials 



(BEDT-TTF) 2 X, with X being I 3 , Cu[N(CN) 2 ]Br or 



2 



Cu(SCN) 2 |33j. Two spin-1/2 kagome-like materi- 
als have been also recently reported, the Volborthite 
Cu 3 V 2 07(OH) 2 2H 2 [34| and the kagome-staircase 
compounds Ni 3 V 2 8 and Co 3 V 2 8 |35| . 

In the present work, we consider the spin-1/2 Heisen- 
berg antiferromagnet on a lattice that interpolates be- 
tween the triangular and the kagome ones. The Hamil- 
tonian is 



H = J^S, 

(ij) 



(if) 



V > 



(1) 



where Sj are spin- half operators, (ij) denotes nearest- 
neighbor bonds belonging to the kagome lattice, and (if) 
are the remaining bonds connecting kagome and non- 
kagome sites. A scheme is indicated in Fig. ^ In this 
way, J' = J corresponds to the usual triangular lattice, 
while J' = defines the kagome one. Our aim is to start 
from the triangular limit and to investigate the stability 
of the ordered state as J' /J decreases. This is sensible 
because the classical Neel state on the triangular lattice is 
compatible with the expected y/3 x \/3 classical ordering 
on the kagome antiferromagnet. The first to study this 
model were Zeng and Elser [23, who performed a spin- 
wave analysis and concluded that, for spin-1/2 particles, 
the ordered state could be stable for J' down to J'/ J ~ 
0.2. Very recently, this model has been investigated with 
a coupled cluster treatment |36| . This technique is based 
on the three-sublattice structure characterizing the 120°- 
Neel order of the triangular lattice, which is found to 
break down very close to J' = 0, indicating the instability 
of the magnetic ordered state very close to the kagome 
limit. 

In this paper, we tackle this problem using variational 
approaches, and exact diagonalization of small clusters. 
In particular, in Section II, we employ the so-called fixed- 
node (FN) technique [HI m order to improve the accu- 
racy of a wave function with long-range antiferromag- 
netic order previously introduced in the pure triangular 
case 0, . This technique, is usually used in the con- 
text of quantum Monte Carlo simulations as a method to 
approximate the Hamiltonians affected by sign-problem 
instabilities |3 1| and obtain exact ground-state properties 
of the corresponding "effective Hamiltonian" , no longer 
affected by the sign problem. Here we use the FN method 
to define a variational state with long-range antiferro- 
magnetic order, and we check its accuracy in describing 
the ground state of the model as J'/ J is reduced through 
a direct comparison with exact diagonalization results on 
the 6x6 cluster. In order to detect any indication of a 
change in the nature of the ground state approaching 
the kagome limit, in Section III we analyze the structure 
of the low-energy spectra as a function of J'/ J, using 
the exact diagonalization of the Hamiltonian on small 
clusters. Section IV is finally devoted to summary and 
conclusions. 




FIG. 1: The depleted triangular lattice. Filled and empty 
circles are the kagome and non- kagome sites; solid and dashes 
lines indicate J and J' bonds of the Hamiltonian Q. The 
letters A,B,C label the three different spin directions oriented 
120° apart of the y/3 x y/3 classical Neel state. 



VARIATIONAL APPROACHES 

A fairly accurate representation of the ground state of 
the spin-1/2 Heisenberg antiferromagnet on the triangu- 
lar lattice can be obtained starting from a 120° Neel or- 
dered state and including Gaussian fluctuations by means 
of a Jastrow factor containing two-spin correlations 



\ip v ) = P exp 



■j)Sfs;)\N) 



(2) 



where Po is the projector onto the S z = subspace, \N) 
is the classical Neel state in the xy plane, 



\N) 



H ex P 



2wi / >— ^ 
"3V ^ 



Sf 



tec 



Si 



(3) 



and |x) is an Ising spin configuration specified by assign- 
ing the value of Sf for each site. On the square lat- 
tice case, the classical Neel state reproduces exactly the 
phases of the ground state of the Heisenberg Hamiltonian 
according to the Marshall theorem. |1 3| On the triangu- 
lar lattice, instead, the exact phases of the ground state 
are unknown, and the classical part of the wave function 
does not reproduces them accurately. However, as 
originally suggested by Huse and Elser, [4J a very accu- 
rate ansatz of the ground-state phases can be obtained 
by including three-spin correlation factors of the form: 



T(x) = exp(z/? J2 ^kS*S*S z k 



(4) 



(i,j,k) 



defined by the coefficients 7^ = 0,±1, appropriately 
chosen so as to preserve the symmetries of the classical 
Neel state, and by the variational parameter (3. In par- 
ticular the sum in Eq. @ runs over all distinct triplets 



3 



of sites i, j, k where both i and k are nearest neighbors of 
j, and i and k are next-nearest neighbors to one another. 
The sign factor 7^ = 7^ — ±1 is invariant under rigid 
translations and rotations in real space by an angle of 
120° of the three-spin cluster k, but changes sign un- 
der rotations by 60°. The resulting wave function reads 
therefore: 



J'/J P 



IT 



Ey/j Eq/J 



= Po^O(x)exp(i^^-i)^j 



with the phase factor given by 



(5) 



(6) 



Since the Hamiltonian is real, a better variational wave 
function is defined by the real part of EqEJ 

The two-body Jastrow potential in j£| contains in prin- 
ciple as many variational parameters as the indepen- 
dent distances on the lattice. However, the same level 
of accuracy can be obtained by optimizing separately 
the nearest-neighbor and next nearest-neighbor distances 
and adopting for the longer-range correlations an expres- 
sion based on the consistency with linear spin-wave the- 
ory: 



v(r) 



N 



with 



l + 2 7fc 
1 - Ik 



(>) 



(8) 



7a,. = 2(cosfc x + 2cosfca;/2cos V3fcy/2) and 7700 is a vari- 
ational parameter. For the anisotropic triangular lattice 
(see Fig^) we have optimized separately the nearest- 
neighbor bonds connecting two kagome sites (rji), and 
a kagome site and a non- kagome site as well as the 
next-nearest-neighbor bonds (772)- As a result, the total 
number of variational parameters for the present vari- 
ational wave function is 5. Their optimal values, and 
the corresponding variational energies, are reported for 
TV = 36 in Tab. [I] for various values of J'/J. 

The accuracy of the present long-range ordered wave 
function in the triangular lattice limit (J'/J = 1) has 
been analyzed in detail in Ref . [T3L IT^j . In this limit the 
wave function is known to provide a qualitatively correct 
representation of the ground-state correlations. 00,0, 
Il5| In order to check the accuracy within a larger range 
of J'/J we have compared several variational properties 
with the exact ground-state values calculated by exact 
diagonalization on the largest cluster presently accessible, 
N = 36. As shown in Fig. El for 0.4 < J'/J < 1 both the 
accuracy on the ground-state energy and the overlap with 
the exact ground state, remain approximatively constant 
and equal to the values in the triangular limit (J'/J = !)• 
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1.00 



10.7687 
11.7299 
12.6911 
13.6522 
15.5730 
17.4873 



■12.0799 
12.7120 
13.5026 
■14.3708 
16.2287 
18.1754 



1.0 0.23 -0.73 -0.73 0.055 1.00 -19.4239 -20.1734 

TABLE I: Variational parameters and variational energies for 
the spin-wave wave function JKJ for different values of the 
ratio J'/J on the N = 36 cluster. The Lanczos exact values 
of the energy are also reported for comparison. 



In addition, the wave function provides an accurate 
representation of the phases of the actual ground state 
up to J'/J ~ 0.2, as it can be checked by measuring the 
average-sign 



s gr# (z)V'o (x)} I ipo (x) 1 2 



(9) 



with \ipo) = J2 X , 4'( x )\ x ) (Fig. |2Jl ■ This remarkable fea- 
ture is due to the presence of the triplet term 10} , al- 
lowing to adjust the phases in a non-trivial way, without 
changing the underlying Neel order. For instance, in the 
triangular limit the average sign (overlap) of the wave 
function is 0.733 (0.562) without the triplet term and 
0.932 (0.779) with it. 

Since the variational wave function (0) reproduces ac- 
curately the phases of the ground state its quality can be 
improved by adopting the FN scheme of Ref. [3l| • This al- 
lows one to obtain a new variational wave function, IV'fn), 
defined as the ground state of the so-called FN effective 
Hamiltonian, whose matrix elements, H xx ,, can be con- 
structed starting from the original Hamiltonian and a 
variational guess on the ground-state phases given, in 
our case, by the wave function JSJ, \%jj v ) = "4>v(x)\x): 




H x >, x < 
H x >, x > 



where H x > >x = ip v (x')H x >^/iJj v (x), and 
V(ar) = Yl R *' 



(10) 



(11) 



Indicating with E^ 1 the ground-state energy of the 
FN Hamiltonian, with E { v n = (ipf n \H\ipf n ) / (ip{ n \ipf n ) and 
E v — (ip v \H\ip v ) / (ip v \i/jv) the energy expectation values 
on IV'fn) and \ip v ), respectively, it is possible to show [3l| 
that the following chain of inequalities holds: 

E v > i?o" > E v n > E , 

where Eg is the ground-state energy of H . Hence, the FN 
procedure is granted to produce a wave function with a 
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FIG. 2: Results on the N = 36 cluster. Upper panel: Relative 
error on the ground-state energy, for the spin-wave (empty 
triangles) and FN (full triangles) wave functions. Stars refer 
to the accuracy of the upper bound on the energy given by the 
lowest eigenvalue of the FN Hamiltonian, Eg 1 . Lower panel: 
average sign of both the spin-wave and FN wave function 
(circles) and their overlap with the exact ground state (same 
symbols as above). Inset: antiferromagnetic order parameter 
The circles are the exact ground-state values, while empty 
and full triangles correspond to the results obtained with spin 
wave and FN wave functions, respectively. Lines are guides 
for the eye. 



better variational energy than \ip v ). In addition, also the 
lowest eigenvalue of the FN Hamiltonian, Eq u , gives an 
upper bound of the ground-state energy better than the 
variational energy E v . This is the quantity which is usu- 
ally considered in the quantum Monte Carlo application 
since it is the most directly accessible. 

The FN Hamiltonian is explicitly defined in such a way 
that the matrix H^ x = sgn[ip v (x')]H x ! f x sgn.[ip v (x)] has 
all negative off-diagonal matrix elements. By the Perron- 
Frobenius theorem, the amplitudes of the ground state 
of this matrix, ip{ n (x), have the same sign for all the 
configurations \x). This implies in turn that the FN 
ground state, IV'fn), and the starting variational wave 
function, \tJj v ), have exactly the same phases. In fact, 
the amplitudes of the ground state of H*? x , ipf n (x), are 
related to the amplitudes of the ground state of the 
transformed matrix H x f x , ^^(x), by the simple relation 
4>fn(z) — ipfn(x)sgn[ip v (x)). For this reason, the FN wave 
function is expected to provide an accurate description 
of the ground state only when it is constructed with a 
good variational ansatz of the ground-state phases. In 
the present case for 0.2 < J' / J < 1. 

In order to thoroughly check the accuracy of the FN 



wave function, we have exactly diagonalized with the 
Lanczos algorithm the FN Hamiltonian for several val- 
ues of J' I J on the N = 36 cluster. This has allowed 
us to calculate not only the FN energies but also the 
overlap of the FN wave function with the exact ground 
state. As shown in Fig. [21 both the FN upper bounds to 
the ground-state energy, and Eq 1 , are sizably more 
accurate than the simple variational estimate E v . In par- 
ticular, the FN wave functions has a much higher overlap 
than \ij) v ), and its accuracy is almost constant and com- 
parable to the one in the triangular limit down to values 
of J'/ J as small as ~ 0.2. 

We have finally compared the variational, FN, and ex- 
act estimates of the antiferromagnetic order parameter 
for a 120° Neel order, 



,t2 



36- 



M 2 



* , (12) 
N(N + 6) v ' 

where M 2 is the sublattice magnetization squared 0. 
Interestingly, by decreasing J'/ J the exact ground-state 
order parameter remains approximatively constant down 
to J'/ J ~ 0.2 thus indicating a possible destabilization of 
the antiferromagnetic order only very close to the kagome 
limit. In addition, though the variational and FN esti- 
mates of the order parameter are approximatively 10% 
higher than the exact one the same degree of agreement 
is observed in all the range 0.2 < J'/ J < 1. For the 
N = 36 cluster, the FN wave function, based on the vari- 
ational wave function (JSJ, provides a good quantitative 
description of the exact ground state for 0.2 < J'/ J < 1. 
In this range, the variational and the exact expectation 
values of antiferromagnetic order parameter remain con- 
stant and equal to their values in the triangular limit. 
Below J'/ J ~ 0.2, the exact value of the order pa- 
rameter begins to decrease and the accuracy of the our 
Neel ordered wave function quickly degrades, indicating 
a change in the ground-state correlations only very close 
to the kagome limit. 

In order to support further the stability of the antifer- 
romagnetic phase for very small values of J'/ J we have 
extended the variational and the FN calculations of the 
order parameter rrv to much larger sizes by using quan- 
tum Monte Carlo techniques. [3?J As shown in Fig. [3] 
even at a value of J 1 / J — 0.2, the lowest coupling ra- 
tio when the variational wave function is expected to be 
accurate (see Fig|2J), the order parameter remains siz- 
ably larger than the triangular J'/ J = 1 case. This indi- 
cates that, as we increase the size at fixed ratio J' /J < 1, 
the stability of the ordered phase, already evident in the 
exact diagonalization in the N = 36 cluster (see Fig|2J, 
becomes more and more clear. We expect that this quali- 
tative behavior, confirmed both by the exact diagonaliza- 
tion, and even more strongly by the variational and the 
FN approaches on larger sizes, is a genuine feature of the 
model, even though the quantitative results that we have 
obtained by quantum Monte Carlo may be affected by 
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TABLE II: Ground-state energy per spin for different values of 
the J' I J ratio obtained with the spin- wave wave function (e„), 
the FN technique (eo 1 ) and the estimated ground-state energy 
(e^ 11 ) obtained by assuming that the FN error in the energy is 
weakly size dependent (see text). The data for J'/J = are 
normalized to give the energy per spin on the kagome lattice. 
Uncertainties are of the order of one unit on the last digit. 
The coupled cluster method results, from Tab. I and Fig. 2 
of Ref. |3q| , are also shown for comparison. 
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a sizable error. This feature may appear rather surpris- 
ing, as the quantum fluctuations should increase for small 
J'/ J and should tend to destabilize the ordered phase, as 
expected for instance within spin- wave theory. 19] How- 
ever, the wave function that we have used is consistent 
with spin wave theory in the large spin limit, and there- 
fore, since also at the variational level the value of m) 
increases, we conclude that the quantum fluctuations are 
not very accurately described within a method that is 
not controlled by the variational principle (the large spin 
limit), at least in the region of small J'/ J. 

A similar size-scaling analysis can be carried out to 
estimate the ground-state energy per spin in the thermo- 
dynamic limit as illustrated in Fig. 0] The extrapolated 
energies are listed in Tab. [D] Since the FN energy error 
is known exactly up to the 6x6 cluster, we can estimate 
the ground-state energy in the thermodynamic limit, by 
adding the 6x6 correction to the infinite-size FN es- 
timate. The corresponding extrapolated values of the 
ground-state energy shown in the Tab^ represent rea- 
sonable benchmark values for this quantity (or at least 
good upper bounds), as the FN error is expected to in- 
crease weakly for larger sizes even for a good variational 
ansatz. This behavior has been verified up to 6 x 6, and 
is also very reasonable to expect in general for an approx- 
imate variational calculation. 

In the triangular limit, J'/J = 1, our estimated 
ground-state energy coincides with the corresponding one 
obtained with the coupled cluster method j3|| by ex- 
trapolating in the size of the clusters (m — > oo, with 
the notations of Ref. H(|). This value is instead slightly 
lower than the extrapolated results based on small clus- 
ters (N < 36) obtained in Ref.0, giving e = -0.5445. 
A good agreement between the coupled cluster for m = 6 
and the FN method is also seen for values of J'/ J down 
to 0.2 (Tab. E3J. Furthermore, our extrapolated ener- 
gies, based on the error for the 6 cluster and the FN 
energy, remain lower than the m = 6 coupled cluster re- 
sult by a similar amount, suggesting that our variational 
approach remains enough accurate also in this region (un- 



FIG. 3: Size scaling of the antiferromagnetic order parameter 
for the spin- wave (empty symbols and dashed lines), and the 
FN (full symbols and dotted line) wave functions: J'/J — 0.2 
(triangles), J'/J = 1.0 (squares). 

fortunately the extrapolations m — > oo are not given in 
Ref.[3(| for < J'/J < 1). Instead in the kagome limit 
our variational ansatz should be clearly less accurate, as 
the ground state is believed to be a spin liquid with no 
magnetic order in the thermodynamic limit, i.e. qualita- 
tively different from our initial variational guess given by 
Eq. J2J). This may explain why in this case our energy 
estimate is slightly higher than the m — > oo coupled clus- 
ter result. We note, however, that our values reported in 
Tab. (flT)) represent in general reasonable upper bounds 
for the energy, as they are obtained by a rigorous varia- 
tional method such as the FN one. 

LOW-ENERGY EXCITATIONS 

An effective method to investigate the possibility of 
magnetic order in spin systems is to analyze the struc- 
ture of the spectrum of finite-size samples, following the 
strategy of Refs. 0,0 EH- The tendency toward antifer- 
romagnetic order in the thermodynamic limit manifests 
itself in finite-size clusters through the fact that the spin 
excitations with the lowest energies can be described by 
the effective Hamiltonian of a "quantum top" . Within 
such a description, the energy of the lowest levels in the 
different subspaces labeled by total spin S, can be ap- 
proximated by, 

E N (S)-E N (0) = ^§±^, (13) 
zi N 

where In is the inertia of the top, which is an extensive 
quantity. Hence, the plots of the lowest energy levels 



6 



-0.31 




FIG. 4: Size scaling of the ground-state energy per spin for 
J' /J = 0.2: Spin- wave wave function (empty triangles), FN 
(full triangles) 



as functions of S(S + 1), have the appearance of a "Pisa 
tower" with a slope that decreases as N increases. An im- 
portant property of the states of the "Pisa tower" (indi- 
cated in Ref. quasi degenerate joint states) is that they 
belong to irreducible representations of the point group 
that are compatible with the symmetry of the ordered 
state. In the case of the -s/3 x V3 order, these are the rep- 
resentations labeled as Ti, T2, T3 of the Cs v group, which 
correspond to [k = 0,7^* = \J/,7£ 2 ,r/3* = *,cr x * = 
*];[k = 0,7^* = -*,^ /3 * - ^o-,.* - tf];[k = 
Q,TZ^ = *,^ 2 7r/3* = *,o"x* = *], respectively. We 
have followed the same notation of Refs. ||| 0, Ha|, 72-0 
denoting a rotation of 4>, cr x , a reflection with respect to 
a mirror plane with the normal pointing along x, and be- 
ing Q = (2tt/3, — 27r/3), the corner of the Brillouin zone 
of the kagome lattice. 

There are not so many finite periodic clusters, which 
can be exactly diagonalized and that interpolate be- 
tween the triangular and kagome lattices without frus- 
trating the antiferromagnetic order. The smallest ones 
are N = 12 and N = 36, which evolve, respectively, 
toward N = 9 and N = 27 at the pure kagome limit 
(</' = 0). The unit cell of the interpolating lattice con- 
tains 2x2 unit cells of the pure triangular lattice, im- 
plying a reduction of the translation operations of the 
periodic cluster of a factor 4 with respect to the number 
of translations at the pure triangular limit. This implies a 
much bigger Hilbert space to be dealt with in the numer- 
ical procedure. In concreteness, even exploiting all the 
available symmetries, the number of states is 42035724 
for N = 36 spins. In this cluster, a complete study of sev- 
eral excited states within different spin sectors for several 
values of J' is prohibitive from the computational point 
of view, while the N — 12 cluster maybe too small. We 
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FIG. 5: Energy spectra as functions of S(S + 1) for a clus- 
ter with N = 12 spins, periodic boundary conditions, and 
J'/J = 0.8, 0.6, 0.4, 0.2 (left to right and top to bottom). The 
crosses (circles) indicate energy levels belonging to symmetry 
representations compatible (incompatible) with the \/3 x \/3 
magnetic order. 
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FIG. 6: Energy spectra as functions of Si for a cluster with 
TV = 16 spins and twisted boundary conditions. Details are 
the same as in Fig. 



have, therefore, included in the analysis the clusters with 
iV = 16 and N = 28 sites (which evolve toward kagome 
clusters with iV = 12 and N = 21, respectively), by in- 
troducin g; tw isted boundary conditions as explained in 
Refs. |t! l26j|. The latter are equivalent to rotate the lo- 
cal frame in (±27r/3, =f27t/3) at each translation along a 
unitary lattice vector. This procedure restores the oth- 
erwise frustrated antiferromagnetic order in the x, y spin 
plane but obviously breaks the spin-rotational symmetry 
as well as some symmetry operations of the point group. 
For this reason, in this law like that expressed 

in (|13fl . with the replacement S(S + 1) — > S%, should 
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FIG. 7: Lowest energy levels within the different S% sub- 
spaces for N — 28 and twisted boundary conditions. Details 
are the same as in Fig. 



be obeyed by the quasi degenerate joint states, as the to- 
tal S z along the z axis, cannot couple with the remaining 
total spin components in the clusters with twisted bound- 
ary conditions. Results are shown in Figs. 03 and for 
N = 12,16, and 28, respectively. Figures and El show 
the full spectra, obtained by diagonalizing all the blocks 
of the Hamiltonian matrix. In the case of N = 28 (Fig. 
|7J), we have not diagonalized the full Hamiltonian but 
obtained the lowest eigenvalues within all the subspaces 
using the Lanczos algorithm. 

We have plotted the states with symmetries compatible 
(non-compatible) with the antiferromagnetic y/3 x y/3 
order with crosses (circles). The first striking feature 
is the fact that only a subset of the states forming the 
"Pisa tower" in the triangular limit remains aligned as 
J'/ J decrease. Those states build up, so to say, a "small 
Pisa tower", whose slope decreases with J'. 

In order to analyze in more detail the behavior of the 
states along the "small Pisa tower" we show zooms of 
the low-energy and low-spin sector of the spectra for the 
N = 12, and N = 28 clusters in Figs. 03 and M re- 
spectively In the case of N = 12, shown in Fig. [H] 
this small set contains levels that belong to representa- 
tions compatible with the \/3 x \/3 order and they remain 
aligned within all the range 0.1 < J'/ J < 1. Signatures 
of departure from that behavior are observed for very 
low J'/ J = 0.05, where the lowest level with 5 = 2 devi- 
ates from the line of the "small Pisa tower" while a state 
with a symmetry not compatible with the antiferromag- 
netic order becomes quasi degenerate with it. In Fig. |S| 
the data for the related N = 9 kagome lattice are also 
shown. In that similar deviation of the state with 

5 = 5/2 from the line connecting the lowest levels in the 
5=1/2 and 5 = 3/2 sectors is observed but within a 
scale which is an order of magnitude larger in the energy 
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FIG. 8: Detail of the low-lying energy levels close to the 
kagome limit for a cluster with N = 12 spins and periodic 
boundary conditions. The lower panel shows the behavior of 
the low energy levels at the kagome limit in a cluster with 
N = 9 sites. Details are the same as in Fig. Qj] 
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FIG. 9: Detail of the low-lying energy levels close to the 
kagome limit for a cluster with N = 28 spins and twisted 
boundary conditions. Details are the same as in Fig. 



axis. Another contrasting feature that comes from the 
comparison between J' = and J' ^ is that in the 
pure kagome limit, there are two states with 5 = 1/2 
(one of them belonging to the manifold of the degenerate 
ground state) whose symmetries are not compatible with 
the -\/3x y/3 order that have energies within the spin gap. 
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FIG. 10: Gap to the lowest excitation with AS Z ~ 1 for TV = 
12 (circles), N = 28 (triangles), and (in the Ti representation) 
TV = 36 (squares). The open circle and square at J ' = 
indicate the value of the spin gap for the TV = 9 and TV = 27 
kagome cluster (corresponding to the TV = 12 and TV = 36 
depleted triangular clusters) . The cross at J' — corresponds 
to a kagome cluster with TV = 36. Inset: The ground-state 
energy as function of J' (open diamonds). The filled diamond 
indicates the ground-state energy of the kagome cluster with 
TV = 27 spins. 



Instead for J'/J = 0.05 the spin gap is clean from such 
states. In addition, the gap to the lowest excitation with 
AS = 1 does not evolve with continuity for J' / J — ► 
to the corresponding gap in the kagome limit (see also 
Fig.EJ. 

Similar remarks apply to the 16- and 28- sites clusters 
shown in Figs. and In these cases, see e.g., Fig. 
the "small Pisa Tower" contains a larger number of states 
but the deviation from a perfect alignment is observed 
already at J'/J < 0.2, where states within subspaces 
not compatible with the V3 X V3 order appear at low 
energies. In particular, for S z = 1 one of those states 
is almost degenerate with the one belonging to the Pisa 
tower while for S z = 2 it is well below it. As in the 
TV = 12 cluster, also for TV = 16 and 28 the gap to 
the first spin excitation monotonically decreases with J' 
(Fig.EJ. 

To complete the analysis, we have computed the 
ground-state energy (Eq) and the lowest eigenvalue in 
the subspace with S = 1 within the representation Ti of 
C^v (Es=i,Ti), in the periodic lattice with TV = 36 sites. 
The latter does not actually correspond to the lowest en- 
ergy S — 1 excitation in the pure triangular limit, which 
is in the subspace corresponding to T2- In any case, the 
energy difference A = £'s = i j r 1 — Eq is, in general, an up- 
per bound for the spin gap of the cluster in all the range 
< J'/J < 1. The behavior of A as a function of J' 



is shown in Fig. ^] As in the clusters previously con- 
sidered, the gap decreases rapidly as J' decreases, and, 
in particular, it does not evolve continuously for J' — * 
towards the value of the corresponding kagome lattice 
(TV = 27, indicated with an open square in Fig. EH - In- 
stead, it tends to a value which is even smaller than the 
magnitude of the spin gap for the TV = 36 kagome lat- 
tice (indicated by a cross in the same figure). This is in 
contrast with the behavior exhibited by the ground-state 
energy, which evolves smoothly toward the value of the 
TV — 27 kagome lattice (see the inset of Fig. II 01) . 

The behavior of the spin gap and the spectra strongly 
suggests the closing of the spin-gap in the thermody- 
namic limit within the whole range of < J'/J < 1. In 
fact, such gap is known to close in the thermodynamic 
limit at the triangular point (J'/J =1) 0,0,0 

and it 

decreases systematically by decreasing J 1 in all the clus- 
ter considered as the low-energy scale for the lowest-spin 
excitations is set by J', i.e., by the states of the "small 
Pisa tower" . This can be understood in a simple un- 
corrected framework where for J' < J the lowest spin- 
excitations are clearly obtained through a spin-flip on a 
non-kagome site, with an energy cost 6J', compared to 
a cost 4 J + 2 J' for a spin flip on a kagome site. As the 
number of non-kagome spins is just a fraction of the to- 
tal number of sites, on a finite cluster such a mechanism 
would apply only to the lowest spin excitations: hence 
the reduction of the number of states belonging to the 
"Pisa tower" with respect to the triangular case. 

This simple picture immediately suggests that the na- 
ture of the spin excitations is intrinsically different for 
finite J' and in the pure kagome limit. In the latter case, 
in fact, the spins on the non-kagome sites do not belong to 
the Hilbert space of the model and such low-energy spin 
excitations are not possible. For this reason by turning 
on the J', the spin gap has a finite discontinuity and it 
is reasonable to expect the model to be unstable under 
the perturbation introduced by J' bonds. This is also 
confirmed by the analysis of the low-energy excitations 
for small values of J'/J (see Fig. |HJ) 

One step further in the analysis of the nature of the 
ground state leads to the question whether or not the lat- 
ter remains ordered down to the kagome limit. The struc- 
ture of the spectra of the small TV = 12 cluster indicates 
that could well be the case, while those of TV = 16 and 
TV = 28 suggests some change in the nature of the ground 
state around J'/J ~ 0.2. It is worth recalling, however, 
that the results for the two latter clusters have been ob- 
tained by using twisted boundary conditions which could 
lead to weaker signatures of ordering. On the other hand, 
in two dimensions, it seems difficult to conceal the pos- 
sibility of the closing of the spin gap with the absence of 
some kind of magnetic order. A simple geometrical anal- 
ysis reveals that the non-kagome spins, which we have 
already identified as the responsible of the lowest-energy 
excitations, form themselves a triangular lattice with a 
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cell parameter twice the size of the one of the usual tri- 
angular lattice. Therefore, a possible scenario for the 
evolution of the ground state as J' decreases could be 
that at some point a crossover takes place from a mag- 
netic order in the usual %/3 x \/3 pattern to a magnetic 
order with a similar pattern but in the triangular lattice 
of the non-kagome spins. Both kinds of order are com- 
mensurate and the ground state could undergo a smooth 
evolution from one to the other. 



SUMMARY AND CONCLUSIONS 

We have investigated the low energy properties of a 
J' — J triangular lattice that interpolates between the 
usual triangular (J' /J = 1) and kagome (J' = 0) lat- 
tices. To this end, we have used a variational approach 
based on a FN wave function accurately describing the 
ground state in the triangular (J' /J = 1) limit, and exact 
diagonalization techniques. 

We have analyzed the quality of the approximation to 
the exact ground state, provided by the FN technique in 
a periodic cluster with N = 36 sites and then extended 
the calculation up to TV = 144 sites by using quantum 
Monte Carlo. We have found that such a wave function, 
describing a state with a v3 x \/3 Neel ordered phase, is 
very close to the exact one for 0.2 < J' /J < 1 and that 
in this range its accuracy is almost independent of the 
J'/ J ratio. Consistently, the ground-state expectation 
value of the antiferromagnetic order parameter remains 
approximately constant down to J'/ J ~ 0.2. Below this 
value the order parameter is suppressed and the quality 
of the FN wave function degrades, suggesting a change in 
the nature of the ground-state. This is also confirmed by 
the analysis of the low-energy spectra on small clusters, 
showing some signatures of instability of Neel ordering 
for J' / J < 0.2. However, at the same time, such analysis 
also indicates quite clearly that the low-energy scale for 
spin excitations is set by J' and that in particular the 
spin-gap has a finite discontinuity at J' / J — 0, the triplet 
excitations being gapless for any non zero value of J'/ J. 
This would imply that the kagome disordered phase is 
unstable against a slight perturbation tending to restore 
the z = 6 coordination number of the triangular lattice. 
A possible scenario is that within the large amount of 
singlets that are quasi degenerate with the ground state 
in the kagome clusters, the one corresponding to the -\/3x 
V3 state is favored by some kick produced by J'. 

Our results are in disagreement with the predictions 
of the spin- wave theory of Ref. 0]j indicating a pro- 
gressive reduction of the antiferromagnetic ordering for 
J'/ J < 1 and a complete melting of the V3 x v^3 order 
for (J'/ J) c ~ 0.2. Instead, our conclusions are closer to 
those of the recent coupled cluster treatment of Farnell et 
al. |36j providing evidence for instability of the antiferro- 
magnetic order very close or possibly at the kagome point 



((J'/J) c = 0.0 ± 0.1), and the existence for small J' / J 
of a regime whose correlations are very different to those 
of the triangular antiferromagnet. Further investigations 
are necessary to clarify the nature of this regime. 

ACKNOWLEDGMENTS 

L. A. thanks Prof. Fulde for his hospitality as well 
as the the Alexander-von-Humboldt Stiftung and CON- 
ICET, Argentina for the support. This work was par- 
tially supported by INFM-PRA MALODI, and by NSF 
under Grant No. DMR02-11166 (L.C.). 



[1] Permanent address: Departamento de Ffsica, Universi- 
dad de Buenos Aires, Ciudad Universitaria Pabellon I, 
(1428) Buenos Aires, Argentina. 

[2] P. W. Anderson, Mater. Res. Bull. 8, 153 (1873). 

[3] P. Fazckas and P. W. Anderson, Philos. Mag. 30. 423 
(1974). 

[4] D. A. Huse and V. Elser, Phys. Rev. Lett. 60, 2531 
(1988). 

[5] R. Singh and D. Huse, Phys. Rev. Lett. 68, 1766 (1992). 
[6] B. Bernu, P. Lecheminant, C. Lhuillier and L. Pierre, 

Phys. Rev. Lett 69 2590 (1992). 
[7] B. Bernu, P. Lecheminant, C. Lhuillier and L. Pierre, 

Phys. Rev. B 50 10048 (1994). 
[8] P. W. Leung and K. J. Runge, Phys. Rev. B 47, 5861 

(1993). 

[9] K. Yang, L. K. Warman, and S. M. Girvin, Phys. Rev. 

Lett. 70, 2641 (1993). 
[10] P. Sindzingre, P. Lecheminant, and C. Lhuillier, Phys. 

Rev. B 50, 3108 (1994). 
[11] T. Jolicoeur, E. Dagotto, E. Gagliano and S. Bacci, Phys. 

Rev. B 42, 4800 (1990). 
[12] P. Azaria, B. Delamonte, and D. Mouhanna, Phys. Rev. 

Lett 70, 2483 (1993). 
[13] L. Capriotti, A. E. Trumper, and S. Sorella, Phys. Rev. 

Lett 82, 3899 (1999); A. E. Trumper, L. Capriotti, and 

S. Sorella, Phys. Rev. B 61, 11529 (2000); L. Capriotti, 

Int. J. Mod. Phys. B 15, 1799 (2001). 
[14] L. Capriotti, PhD. Thesis, SISSA, Trieste, Italy (2000). 
[15] M. Boninsegni, Phys. Rev. B 52, 15304 (1995). 
[16] D. A . Huse and A. D. Rutenberg, Phys. Rev. B 45, 7536 

(1992) . 

[17] J. T. Chalker, P. C. W. Holdsworth, and E. F. Shender, 

Phys. Rev. Lett. 68, 855 (1992). 
[18] J. N. Reimers and A. J. Berlinsky, Phys. Rev. B 48, 9539 

(1993) . 

[19] C. Zeng and V. Elser, Phys. Rev. B 42, 8436 (1990). 
[20] J. T. Chalker and J. F. G. Eastmond, Phys. Rev. B 46, 
14 201 (1992). 

[21] P. W. Leung and V. Elser, Phys. Rev. B 47, 5459 (1993). 
[22] N. Elstner and A. P. Young, Phys. Rev. B 50, 6871 

(1994) . 

[23] C. Zeng and V. Elser, Phys. Rev. B 51, 8318 (1995). 
[24] A. Chubukov, Phys. Rev. Lett. 69, 832 (1992). 
[25] L. O. Manuel, A. E. Trumper, C. J. Gazza, and H. A. 
Ceccatto, Phys. Rev. B 50, 1313 (1994). 



10 



[26] P. Lecheminant, B. Bernu, C. Lhuillier, L. Pierre and P. 

Sindzingre, Phys. Rev. B 56 2521 (1997). 
[27] C. Waldtmann, H.-U. Everts, B. Bernu, C. Lhuillier, P. 

Sindzingre, P. Lecheminant, L. Pierre, Eur. Phys. Journal 

B, 2 501 (1998). 
[28] M. Mambrini and F. Mila, Eur. Phys. J. B 17 651 (2000). 
[29] C. Lhuillier and P. Sindzingre, "Quantum Properties 

of Low-Dimensional Antiferromagnets" , Editors: Yoshi- 

tami Ajiro and Jean-Paul Boucher, Kyushu University 

Press, 111 (2002). 
[30] G. Misguich and C. Lhuillier, e. print cond-mat/0310405 
[31] D. F. B. ten Haaf, H. J. M. van Bemmel, D. F. B. ten 

Haaf, J. M.J. van Leeuwen, W. van Saarloos, and D. M. 

Ceperley, Phys. Rev. B 51, 13039 (1995). S. Sorella and 



L. Capriotti, Phys. Rev. B 61, 2599 (2000). 
[32] K. Takada et al, Nature 422, 53 (2003). 
[33] R. H. McKenzie, Science 298 820 (1997). 
[34] Z. Hiroi, M. Hanawa, N. Kobayashi, M. Nohara, H. Tak- 

agi, Y. Kato and M. Takigawa, J. Phys. Soc. Jpn. 70 

3377 (2001). 

[35] N. Rogado, G. Lawes, D. A. Huse, A. P. Ramirez, R. J. 

Cava, Solid State Commun. 124, 229 (2002). 
[36] D. J. J. Farnell, R. F. Bishop, and K. A. Gernoth, Phys. 

Rev. B 63, 220402 (2001). 
[37] M. Calandra and S. Sorella Phys. Rev. B. 57, 11446 

(1998). 



